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The calculation of one loop integrals at finite temperature requires the evaluation of 
certain series, which converge very slowly or can even be divergent. Here we review a new 
method, recently devised by the author, for obtaining accelerated analytical expressions for 
these series. The fundamental properties of the new series are studied and an application to 
a physical example is considered. The relevance of the method to other physical problems 
is also discussed. 
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1 Introduction 

In two recents papers [DUI by the author the problem of obtaining fastly conver- 
gent series starting from slowly convergent ones was considered: in pp new series 
representations for the Riemann, Hurwitz and more general zeta functions were 
discovered; in p! some of the results of were later applied to the calculation of 
one-loop integrals at finite temperature. 

The method that we have devised in is based on the Linear Delta Expan- 
sion (LDE), a non-perturbative technique which has been applied with success in 
the past to a large class of problems (see |3] and references therein). The idea is to 
interpolate a given problem, which cannot be solved analytically, with a solvable 
problem, depending upon one or more arbitrary parameters. By then identifying a 
"perturbation" and performing an expansion in such term, a series depending upon 
the arbitrary parameters is found. Finally, by applying the "Principle of Minimal 
Sensitivity" (PMS)[4| the optimal values for the arbitrary parameters can be ob- 
tained. Notice that the expansions that are obtained following such procedure are 
really non-perturbative, since they do not involve any parameter in the problem 
to be small. 

This paper is organized as follows: in Section |5] we describe the method and 
apply it to obtain a new series representation for a generalized zeta function; in 
Section O we consider the application of the method to the calculation of one-loop 
integrals at finite temperature; finally, in Section^ we draw our conclusions. 
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2 The method 

We consider the generalized zeta function defined as 

wliere us > 1. Special cases of cq. are the Riemann zeta function, which is 
obtained for u = I and ^ = 1 and the Hurwitz zeta function which corresponds to 
u = l. 

We will now prove the following 
2.1 Theorem 

It is possible to write ({u, s, ^) as 

T{k + s) 

k=0 

where 



and > > 0). 

2.1.1 Proof 

Our proof is based on the identity 



1 °° 1 1 1 



e ^i"'" + (1 + A(?i)) 
where 

AW ^ (5) 

Provided that A^ > and ^ > —1, it is possible to fulfill the condition 
|A(n)| < 1 and therefore apply the binomial theorem to obtain 



1 _^r{k + s) 



k 



(1 + A{n)y ns) k\ 

Using this result in eq. we have 

7r - 1 >^ ^ r(fc + ^ ( k\ A^(^-^) {~^y 

^^^'^^ - + 1^ 1^ Y{s) k\ ^ \ j J {I + X'Y+^ n<^+i) - ^ ' 

^ n=l k=0 ^ ' j=0 ^ ■' ^ ^ ' 
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By performing the sum over n one obtains 

which completes our proof. 

Having proved our fundamental result, eq. lj2Jl, we now discuss some of the 
properties of the new series: we first stress that the series that we have obtained 
is independent of the arbitrary parameter A although A appears explicitly in the 
expression. This happens because eq. ^ is independent of A and it has just been 
proved that our new series, eq. (0) converges to eq. Q provided that A^ > (^ — 1)/2 
and ^ > — 1. In other words we can say that eq. |(5J describes a family of series, 
each corresponding to a different value of A and each converging to the same series, 
eq. (P). 

We can now consider the "truncated" series 

C(^)(A,«,5,0 = -^ + J2^%7^'^k{\u,s,M') 
^ k=o ^ ^-"^^ 

obtained by restricting the infinite sum to the first iV + 1 terms. Obviously 
^(•^)(A, u, s, ^) depends upon A as a result of having neglected an infinite num- 
ber of terms. We can use this feature to our advantage and fix A so that the 
convergence rate of the series is maximal. 

The proper value of A is chosen using the PMS 4 , which amounts to find A 
fulfilling the equation 

^Cw(s,A)=0. (7) 

A straightforward mathematical interpretation of this condition is that the value 
of A complying with this equation also minimizes the difference ^ 

- — 1 ^ 

C(m, S,,^) - C(Ar)(A,M, S,,^)J . (8) 

To lowest order, which corresponds to choosing iV = 1, one obtains the optimal 
value 



^PMS-\I^ > (9) 



which can be used as long as X%ms > ^ Notice that, since XpMs depends 

upon ^ then C{n){'^pms , u, s, ^) will not be a polynomial in ^: on the other hand, if 
we had chosen A = 0, then C(Ar)(Oj '"i s, ^) would be a polynomial of A^*'* order in ^; 
in that case however the convergence of the series would be strictly limited to the 
region ^ < 1. For this reason we will refer to our accelerated series corresponding to 
ApA/s and to A = as being "nonperturbative" and "perturbative" respectively ^). 

^ ) Typically the labels "perturbative" and "nonperturbative" are used to refer to expansions in 
some coupling constant. As we will see shortly in a physical application the parameter ^ can be 
related to an inverse temperature and not to a coupling constant. 
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Fig. 1. ((N) (-^1 2, 1, 1) as a function of A for A*" = 1, 3, 5. The dotted line corresponds to 
the exact value C(2, 1, 1) = (1 + 7rcoth7r)/2. 



As an illustration wc have compared in Fig. ^tlie exact result for C(2, 1, 1) = 
(1 + 7rcotli7r)/2 (the dotted line) with the approximation C(Ar)('^i 2, 1, 1), with N — 
1,3,5. The maximum of the blue curve corresponds to the value of A given by 
eq. ©. 

In Fig.|21we have plotted the difference |C(2, 1, 1) — J2n=o (^I^TT) I ^ ^ function 
of N and compared them with the difference |C(2,1,1) — (■(^•)(Ap7\/5, 2, 1, 1)| for 
= 1, 5, 10. The first 10 terms of the accelerated series, obtained by using eq. (O 
with A given by eq. © are able to provide the same accuracy of the first 10^ terms 
of the original series, eq. (Q). 

3 Application: one loop integrals at finite temperature 

One-loop integrals at finite temperature have the typical form 

J{m,a,b) = T^i'^yf^^ J- , (10) 



where fi is an energy scale brought in by dimensional regularization, a and b are 
integers and Z? = 3 — 2e is the number of spatial dimensions. w„ = 2'KnT are 
the Matsubara frequencies which appear in the imaginary time formalism. We will 
follow the notation set in jH] and write the expression for a general one-loop integral 
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Fig. 2. The difFerence |C(2, 1, 1) — '^^-o (n^+i) I ^® ^ function of A^. The horizontal hnes 
correspond to iC(2, 1, 1) - C(jv) (^^'ms, 2, 1, 1)| for = 1,5, 10. 



at finite temperature in the form 



(fc2) 



2^a 



n = — CX) 



J{m,a,b) = T^^^ 



/Li being the scale brought in by dimensional regularization, a and b being integers 
(a > and 6 > 0). Following [2 El we also define 



\2tiT 



and obtain 



J(Af,a,6) = T (27rT)^+^'^-^'' 217^^ ^^^/^+"^ 5(M, (11) 



(27r)^ T{DI2) T{h) 



where 



n— 1 ^ ' 



(12) 



and l = b — a — D/2 and D — 3 — 2 e. Depending upon the value of I, the series of 
eq. ()12fl could be divergent and therefore need regularization. 
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Eq. H12|l is clearly of the form considered in eq. (jSJ and one can write 

S{M,l) = --1^ + C(2,/,M^) , (13) 



where 



C(2,;,M2) = _l^ + g£^vl.,(A,2,;,M2) (14) 

fc=0 ^ ' 



and 



Once more we stress that although the series all converge to the same result in- 
dependently of A (provided that A > ^^), the partial sums, obtained by truncating 
the series to a finite order will necessarily display a dependence on the parameter. 
Such dependence, which is an artifact of working to a finite order, will also make the 
rate of convergence of the different elements of the family A-dependent. Since it is 
desirable to obtain the most precise results with the least effort, one will select the 
optimally convergent series by fixing A through the "principle of minimal sensitiv- 
ity" (PMS) '7'. The solutions obtained by applying this simple criterion display in 
general the highest convergence rate and, once plugged back in the original series, 
provide a non-polynomial expression in the natural parameters. 

In the present case, the "natural" parameter in eq. Q is Af ^ and to first-order 
the PMS yields 



- M ^iM, 

On the other hand, if the value A = is chosen, then one obtains the series 

C(2,^,M^) . ^ + f:^T^C(2(/ + fc)), (16) 

which, to a finite order yields a polynomial in M. Notice that such series corre- 
sponds to the expansion used in ;6 and converges only for < 1: on the other 
hand the series corresponding to )^p\.jg converges for all values of M. 

In 12] we have studied both series, corresponding to choosing Ap|,jg and A = 
0, showing that the first few terms of the PMS series provide a very accurate 
approximation to the exact result. For example, C(2, 3/2,1) calculated with the 
first 10^ terms in eq. (pj gives the result 1.5124349215 0: the same precision is 
reached with our PMS series with only 27 terms. 



6 



Non perturbative regularisation of one loop integrals at finite temperature 



In |2] we have used this series to evaluate the one-loop self-energy 



I(rn) 



+ 00 

E 



n— — 00 

'Di 



(2^ 



d^k 1 
1 



M2 



J(A/,0,1) 



(17) 



where 



J(M, 0,1) 



-2-2e 



^-3/2-3e ^2e r(-l/2 + e) 



-]i7^+C(2,-l/2 + .,M2) 



(18) 



and d^fc = drj^-i fc^-^ dfc and = 3 - 2e. 

By using our eq. H16|l and retaining only the divergent terms and those inde- 
pendent of e we obtain 

■^2 _2 n m\{3) 



mT T 7X1" 



1 , 47rT2"l 

- + 7 - log 2" 



8(27r)4r2 



5to\(7) 



7mi"C(9) 



1024 7r6r4 ' 32768 7r8 re 262144 TfiOTS ' 

which reproduces the results of [S], which, however, were considered only up to 
order . As anticipated, the formula obtained is a polynomial in m. 
Using the PMS series to first order we obtain 



JpMsiM, 0,1) 



-Kl 



v/C(3) + M2C(5) ' 
C(5)' 



K2 



(20) 



where ki^2 are the constants of integration independent of M: 



«2 = TTT 



1 

12 



- + 7 - log 
e 



C(3)^ 



/47rT2 



2C(3)' 
C(5) 



C(5)^ 



As mentioned before the PMS result is non-polynomial in M and at a given 
order contains more information than the corresponding result for A = obtained 
at the same order. This can be verified by Taylor expanding eq. ()20|l in M and 
comparing it with eq. H19() . which is calculated to order M^: 



mT 

^(^^'0'^) " -17+12^16^ 
5m\(5)^ 



1 , AttT^ 

- + 7 - log — — 

7miOC(5)^ 



32768^8^6^(3) 262144 T8C(3)' 
7 



8(27r)4r2 1024 TT^T^ 
(21) 
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Not only eq. (|2()(l reproduces correctly the terms going as and of the 
perturbativc series but it also provides all the coefficients of the higher order terms 
apart from factors depending on the Riemann function evaluated at odd integer 
values. 

These properties has been used in 1^ to obtain a variational acceleration of the 
series for the thermodynamic potential, which has then been compared with the 
high-temperature series of Haber and Weldon |Sj , showing that the PMS result to 
a finite order provides an accurate approximation valid also to low temperature. 

4 Conclusions 

We have summarized some of the results recently obtained by the author and 
contained in 2 : a given slowly convergent can be transformed into a fastly con- 
vergent series by introducing into the series an arbitrary parameter and by then 
carefully identifying a term in which to expand. This procedure in general provides 
series which converge very fast, typically geometrically, and for all values of the 
parameters. In the present case the natural parameter in the problem is the inverse 
temperature: we have compared our results with similar results in the literature, 
amounting to a high temperature expansion, showing that the present approach 
provides quite precise result at high and low temperature. Possible further appli- 
cations of the ideas exposed in this paper and in (2, include the calculation of the 
Casimir effect, which is presently being carried out [H]. 

The author acknowledges support of CONACYT grant 40633. 
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